Hot Deformation Characteristics and Dynamic Recrystallization Mechanisms of a Semi-Solid Forged AZ91D Magnesium Alloy

Magnesium alloys show great promise in high-speed transport, aerospace, and military technology; however, their widespread adoption encounters challenges attributed to limitations such as poor plasticity and strength. This study examines the high-temperature deformation of semi-solid forged AZ91D magnesium alloy through a combination of experiments and simulations, with a focus on comprehending the influence of deformation conditions on dynamic recrystallization (DRX). The findings disclose that conspicuous signs of DRX manifest in the yield stress curve as strain increases. Additionally, decreasing the strain rate and temperature correlates with a reduction in both yield stress and peak strain, and the activation energy is 156.814 kJ/mol, while the critical strain and peak strain remain relatively consistent (εc=0.66208εp). Microstructural changes during high-temperature deformation and the onset of DRX are thoroughly examined through experimental methods. Moreover, a critical strain model for DRX and a predictive model for the volume fraction of DRX were formulated. These equations and models, validated through a combination of experiments and simulations, serve as invaluable tools for predicting the mechanical behavior and microstructural evolution, which also establishes a foundation for accurately predicting the deformation behavior of this alloy. By analyzing the hot deformation characteristics and dynamic compression mechanism of the newly developed semi-solid forging AZ91D magnesium alloy, a numerical simulation model can be effectively established. This model objectively reflects the changes and distributions of stress, strain, and rheological velocity, providing a scientific basis for selecting subsequent plastic deformation process parameters and designing mold structures.


Introduction
Magnesium and its alloys are widely utilized in industries such as automotive, aircraft, electronics, energy, and biomedical sectors due to their impressive attributes.These include excellent damping properties, specific stiffness, specific strength, specific rigidity, and ease of processing [1][2][3][4].However, challenges arise from the hexagonal close-packed structure, limited room-temperature plasticity, and constraints in processing techniques, corrosion resistance, and cost [5][6][7].Thermal processing above the recrystallization temperature is required in magnesium alloy production, influenced by factors like temperature, strain rate, and deformation, leading to a complex recrystallization process.Scholars focus on the hot deformation behavior and constitutive equations, like the Johnson-Cook model [8], Khan-Huang-Liang model [9], Fields-Backofen model [10], Voce-Kocks model [11], and Arrhenius model [12].The hyperbolic sinusoidal form of the Arrhenius equation, widely used in materials research, is introduced by Sellars and Tagert [12].Duan's study investigates the flow behavior and microstructure of AZ80A magnesium alloy, exploring the influence of deformation temperature and strain rate on dynamic recrystallization (DRX) structure evolution and strain rate sensitivity index [13].Luan conducts isothermal compression tests on AZ31 magnesium alloy to elucidate its deformation characteristics, indicating dynamic flow softening with accurate constitutive model predictions matching observed values [14].Chen analyzes the microstructure and mechanical properties of Mg-Y-Zn alloy, developing processing charts for optimizing hot working and structural control, with predicted outcomes closely matching actual results [15].Additionally, DRX is a vital mechanism for grain refinement and softening during deformation, such as forging, extrusion, stamping, and rolling.Research on magnesium and its alloys has identified various DRX mechanisms, including low-temperature dynamic recrystallization (LTDRX), discontinuous dynamic recrystallization (DDRX), continuous dynamic recrystallization (CDRX), twinned dynamic recrystallization (TDRX), and rotational DRX (RRX).Galiyev's study revealed LTDRX in Mg-5.8Zn-0.65Zralloys, showing a significant grain size reduction and increased hardness at 423 K and a strain rate of 2.8 × 10 −3 s −1 [16].Tan explored CDRX nucleation mechanisms in Mg-9Al-Zn alloy after high-temperature extrusion and torsion experiments [17].DDRX research highlighted distinct nucleation and grain growth stages [18].The research indicates that deformation temperature and strain rate largely determine the sufficiency of dynamic recrystallization in magnesium alloys.By selecting appropriate deformation temperatures and strain rates, different types of DRX can be induced, which eliminates stress concentration caused by work hardening and reduces the force required during the processing of magnesium alloys.However, investigations into DRX in magnesium alloys are limited, warranting further exploration.Establishing a correlation between the material's constitutive model and DRX mechanisms is crucial for understanding deformation behavior.
It is well known that proper setting of process parameters such as extrusion, rolling, and forging can induce DRX and grain refinement in magnesium alloys, thus enhancing their microstructure and mechanical properties [19].However, these conventional plastic deformation methods have several shortcomings.For instance, rolling alloy plates tend to form a strong basal plane texture, which hinders the secondary processing of deformed magnesium alloys and limits their applications [20].The extrusion process tends to form a ring-like basal plane texture, leading to significant anisotropy in the alloy [21].Additionally, the internal DRX is often incomplete, resulting in coarse uncrystallized regions and an uneven grain size distribution.This affects the material's overall performance and results in high waste and a low finished product rate.Although it has been reported that conventional extrusion processes can refine the grain, the extent of grain refinement is generally limited.The development of severe plastic deformation (SPD) technology can further refine the grain size of the material, producing ultra-fine-grained materials [22].However, SPD is highly complex, involving more processing stages and steps.Moreover, SPD technology imposes strict size limitations on the processed materials, restricting its application to smaller workpieces, primarily in laboratory settings, thus preventing large-scale industrial applications.To address the shortcomings of conventional machining processes and SPD technology, a new process combining pre-deformation with extrusion has been proposed.Jiang et al. [23] studied the Mg-1.58Zn-0.52Gd(wt.%) alloy by applying hot pre-forging at 350 • C and 400 • C followed by slow extrusion.They found that the extruded alloys without pre-forging exhibited a bimodal microstructure consisting of fine DRXed grains and coarse unrecrystallized grains with a strong basal texture.After pre-forging, the extruded alloys exhibited a finer, fully DRXed microstructure and a weaker extruded texture.Compared to the extruded non-wrought alloys, the pre-wrought alloys exhibited an excellent combination of properties in terms of strength, ductility, and yield asymmetry.This predeformation process generates a large number of twins, subgrains, and other defect structures in magnesium alloy ingots.These defects provide additional nucleation sites for DRX during subsequent extrusion, promoting recrystallization nucleation and grain refinement, thereby improving the alloy's strength and plasticity.
Pre-deformation combined with subsequent plastic deformation is simpler than the intense plastic deformation method and can produce finer grains than conventional extrusion, improving the alloy's tensile asymmetry [24].However, after pre-forging, the ingot must be reheated.During reheating, the ingot's internal twins and other defects can diminish, leading to grain growth and reducing the effectiveness of the pre-deformation.To maximize the benefits of pre-deformation, it is necessary to explore the metal forming processes following pre-deformation, particularly semi-solid die forging.This process can eliminate the need for a second preheating between pre-forging and extrusion, and it can significantly enhance the molding capability of the semi-solid billet.Considering the high concentration of the second phase within the alloy and its complex structure [25], it is imperative to investigate the thermal deformation and DRX behavior of semi-solid forged AZ91D magnesium alloy to construct a comprehensive material model, encompassing dynamic recovery and a recrystallization rate model tailored to DRX conditions.By analyzing the hot deformation characteristics and dynamic compression mechanism of the newly developed semi-solid forging AZ91D magnesium alloy, a numerical simulation model can be effectively established.This model objectively reflects the changes and distributions of stress, strain, and rheological velocity, providing a scientific basis for selecting subsequent plastic deformation process parameters and designing mold structures.Therefore, it is paramount to research and develop novel magnesium alloy forming methods, investigate the organizational evolution of fabricated parts, establish an intrinsic model, and lay a solid foundation for subsequent simulation research, thereby advancing the application of magnesium alloys across diverse fields.

Experimental Material
The AZ91D magnesium alloy, detailed in Table 1, was produced through semi-solid forging.The melting equipment used for metal casting in this experiment is a custom-built well-type resistance melting furnace (branch torch furance, Luoyang, China).The heating power is 15 kW, the rated temperature is 1000 • C, and the heating rate is 20 • C/min.Additionally, the furnace cover has an open mouth, facilitating easy transport and casting of molten metal.Furthermore, the furnace cover is equipped with a gas inlet to maintain a protective gas atmosphere throughout the melting process.Specifically, the furnace is filled with a mixture of CO 2 and SF 6 gas in a volume fraction ratio of 99:1.During the melting process, a graphite crucible with a 1 kg capacity is used for transfer metal melt.The casting mold, made of H13 steel with a preheating temperature of 300 • C, is coated with water-based graphite as a release agent and has a diameter of 60 mm with a thickness of 30 mm.The melting temperature is set to 750 • C.After the alloy in the crucible is fully melted, the temperature is maintained for 15 min before cooling.As illustrated in Figure 1a, the melt is allowed to cool to 595 • C before being quickly poured.After pouring, the melt is immediately forged under 170 MPa pressure for 3 min.Finally, after the casting is cooled, it is removed and a sample with a diameter of 8 mm and a height of 12 mm is cut from the mid-height of the casting for compression testing.Using a Gleeble-3800 thermal simulation tester (Poestenkill, NY, USA), the hot compression experiments were conducted with the specimen size of Φ6 mm × 9 mm at temperatures of 648 K, 673 K, 698 K, and 723 K, with compression rates of 0.001 s −1 , 0.01 s −1 , 0.1 s −1 , and 1s −1 .Post-compression, specimens (Figure 1c,d) were sectioned along their central axis, meticulously polished, and treated with a 4% nitric acid solution for corrosion.Transmission electron microscopes (TEM, JEOL JEM-2100F, Tokyo, Japan) were also used to characterize the microstructures in this work.

Finite Element Models
As AZ91D magnesium alloy was not available in the Deform-3D material library, the derived equations were added to create the AZ91D material profile.Physical parameters from data sources [12], detailed in Table 2 and measured in SI units, were integrated.For the finite element model in hot compression, the workpiece was considered plastic, and the upper and lower indenters were treated as rigid.Meshing was applied only to the workpiece with a grid resolution of 30,000 elements.To simulate isothermal compression, the workpiece, indenter, and surroundings were adjusted to match experimental conditions.

Stress-Strain Curve
Figure 2 displays the stress-strain curve of AZ91D magnesium alloy, exhibiting a characteristic DRX curve.The curve indicates an increase in flow stress with higher strain rates and lower deformation temperatures.As deformation continues, the flow stress rises significantly, reaching a peak value.Subsequently, it decreases beyond the peak, stabilizing at higher temperatures or lower strain rates.This highlights the heightened sensitivity of AZ91D magnesium alloy's flow stress to variations in strain rate and deformation temperature during hot compression.

Finite Element Models
As AZ91D magnesium alloy was not available in the Deform-3D material library, the derived equations were added to create the AZ91D material profile.Physical parameters from data sources [12], detailed in Table 2 and measured in SI units, were integrated.For the finite element model in hot compression, the workpiece was considered plastic, and the upper and lower indenters were treated as rigid.Meshing was applied only to the workpiece with a grid resolution of 30,000 elements.To simulate isothermal compression, the workpiece, indenter, and surroundings were adjusted to match experimental conditions.

Stress-Strain Curve
Figure 2 displays the stress-strain curve of AZ91D magnesium alloy, exhibiting a characteristic DRX curve.The curve indicates an increase in flow stress with higher strain rates and lower deformation temperatures.As deformation continues, the flow stress rises significantly, reaching a peak value.Subsequently, it decreases beyond the peak, stabilizing at higher temperatures or lower strain rates.This highlights the heightened sensitivity of AZ91D magnesium alloy's flow stress to variations in strain rate and deformation temperature during hot compression.

Constitutive Equation
Throughout the deformation, where the interplay between work hardening and softening occurs, the Garofalo creep equation emerges as a precise descriptor for the relationship between flow stress (σ), strain rate ( .ε), and temperature (T) [26][27][28]: where A 1 , A 2 , A, α, β, n, n 1 are the material parameters; α for the stress level parameter, α = β/n 1 ; A is the structure factor; n is the stress index; Q is the activation energy; R is the standard molar gas constant (8.314J/mol/K).The logarithmic treatment of Equation (1) yields: Per McQueen's findings [29], the peak stress (σ p ) is predominantly chosen for establishing their constitutive equations.As depicted in Figure 3, there is a distinct linear correlation between ln .ε and lnσ, ln Combining Figure 3  Z-parameter is a convenient way of describing the effect of strain temperature and strain rate during the forming process [12,[30][31][32].
Bringing Equation (5) into Equation (1) and taking its logarithms, one can obtain Equation ( 6): ( ) Z-parameter is a convenient way of describing the effect of strain temperature and strain rate during the forming process [12,[30][31][32].
Bringing Equation (5) into Equation ( 1) and taking its logarithms, one can obtain Equation (6): Figure 4 demonstrates a strong linear relationship between lnZ and ln[sinh(ασ)], the n-value (3.60349) differing by only 0.01288 from the computed (3.61637); this minimal discrepancy corresponds to an error rate of no more than 0.36%, affirming the precision and reliability of using the Z-parameter as a reflection of its σ.
According to Equation ( 6), the σ can be expressed as Equation ( 7): Bringing the data into Equation (8) yields the improved Arrhenius equation: Figure 4 demonstrates a strong linear relationship between lnZ and ln[sinh(ασ)], the n-value (3.60349) differing by only 0.01288 from the computed (3.61637); this minimal discrepancy corresponds to an error rate of no more than 0.36%, affirming the precision and reliability of using the Z-parameter as a reflection of its σ.According to Equation ( 6), the σ can be expressed as Equation ( 7): Bringing the data into Equation ( 8) yields the improved Arrhenius equation: Which: ( )

Critical Condition Model for DRX
DRX involves the formation and growth of recrystallized grains in metallic materials, transpiring under favorable deformation conditions-comprising temperature and strain rate.Hence, meticulous control over these parameters proves crucial in attaining the targeted grain size and its associated properties.Therefore, it is necessary to understand the critical condition model for DRX.

Critical Condition Model for DRX
Per the Sellars model, the critical condition model for DRX integrates both peak and critical strain models [33,34].Model expression for peak strain in the Sellars model: where εp is the peak strain, d0 is the initial grain size.Taking Equation (11) into Equation (10), the peak strain model can also be expressed as: ( ) where Q1 is the deformation activation energy; a1, n1, m1, C1 is a material-related parameter.

Critical Condition Model for DRX
DRX involves the formation and growth of recrystallized grains in metallic materials, transpiring under favorable deformation conditions-comprising temperature and strain rate.Hence, meticulous control over these parameters proves crucial in attaining the targeted grain size and its associated properties.Therefore, it is necessary to understand the critical condition model for DRX.

Critical Condition Model for DRX
Per the Sellars model, the critical condition model for DRX integrates both peak and critical strain models [33,34].Model expression for peak strain in the Sellars model: where ε p is the peak strain, d 0 is the initial grain size.Taking Equation (11) into Equation (10), the peak strain model can also be expressed as: where Q 1 is the deformation activation energy; a 1 , n 1 , m 1 , C 1 is a material-related parameter.Since the internal grain sizes of AZ91D magnesium alloy are 57.8 µm, so the effect of d 0 can be neglected.And Equation ( 12) can be expressed as Equation ( 13): From Equation (13), Figure 5a,b, Equation ( 14) can be obtained: From Figure 5, one can get m 1 = 0.2662, Q 1 = 20,514.87814J/mol, and a 1 = 0.004569, so the peak strain model of AZ91D magnesium alloy can be obtained: ( ) From Figure 5, one can get m1 = 0.2662, Q1 = 20,514.87814J/mol, and a1 = 0.004569, so the peak strain model of AZ91D magnesium alloy can be obtained:

Critical Strain Model
Jonas' findings propose that the onset of DRX modifies the internal thermodynamic system of the alloy [35].In the results, an inflection point appears in the θ-σ curve, indicating the critical conditions for DRX initiation.The expression for the work-hardening rate is outlined below:

Critical Strain Model
Jonas' findings propose that the onset of DRX modifies the internal thermodynamic system of the alloy [35].In the results, an inflection point appears in the θ-σ curve, indicating the critical conditions for DRX initiation.The expression for the work-hardening rate is outlined below:  Jonas and Najafizadeh confirmed that the vicinity of the inflection point of the θ-σ curve can be characterized by a third-order polynomial [35]: where A1, A2, A3, B are fitting constants.So: From Equation ( 18), one can see that −(∂θ/∂σ) has a quadratic relationship, and combined with Figure 6, the −(dθ/dσ) − σ curve is made, which is shown in Figure 7. Jonas and Najafizadeh confirmed that the vicinity of the inflection point of the θ-σ curve can be characterized by a third-order polynomial [35]: where A 1 , A 2 , A 3 , B are fitting constants.So: From Equation ( 18), one can see that −(∂θ/∂σ) has a quadratic relationship, and combined with Figure 6, the −(dθ/dσ) − σ curve is made, which is shown in Figure 7.
where A1, A2, A3, B are fitting constants.So: From Equation (18), one can see that −(∂θ/∂σ) has a quadratic relationship, and combined with Figure 6, the −(dθ/dσ) − σ curve is made, which is shown in Figure 7.According to Figure 7, the σ c , σ p , ε c , and ε p for different deformation conditions can be obtained, and the σ c − σ p linear fit and ε c − ε p linear fit are plotted in Figure 8.It is evident that the two fits align perfectly, aligning with the Sellars model's principle that the ratio of ε c to ε p should fall within the range of 0.6 to 0.95.In this study, the relationship equation between σ c and σ p , ε c and ε p can be expressed as Equation (19).According to Figure 7, the σc, σp, εc, and εp for different deformation conditions can be obtained, and the σc − σp linear fit and εc − εp linear fit are plotted in Figure 8.It is evident that the two fits align perfectly, aligning with the Sellars model's principle that the ratio of εc to εp should fall within the range of 0.6 to 0.95.In this study, the relationship equation between σc and σp, εc and εp can be expressed as Equation (19).

DRX Kinetic Model
Currently, the widely employed DRX kinetic model, based on the Avrami equation, remains a prominent choice.In numerous hot simulation experiments, time is introduced as a variable to formulate the following DRX kinetic model [36]: Yada modified DRX's kinetic model by replacing time with strain [37].

DRX Kinetic Model
Currently, the widely employed DRX kinetic model, based on the Avrami equation, remains a prominent choice.In numerous hot simulation experiments, time is introduced as a variable to formulate the following DRX kinetic model [36]: Yada modified DRX's kinetic model by replacing time with strain [37].
where X drx is the percentage of DRX, x and z are material-related parameters, and f (ε) is a strain-related quantity.Furthermore, the DRX kinetics model can be expressed as: where X drx is the volume fraction of DRX, Q 5 is the thermodynamic activation energy of DRX, .
ε is the strain rate, ε 0.5 is the strain corresponding to 50% of DRX, and β d , a 10 , k d , h 5 , m 5 , and c 5 are material correlation coefficients.
Taking logarithms of Equation ( 22) and writing a 10 ε p as ε c , one can obtain Equation ( 24): Figure 9a illustrates that DRX occurs after the material reaches σ p , and the curve descends to the steady-state stress σ ss , indicating material softening.Here, the degree of stress reduction from σ p is denoted as σ p − σ exp , and the maximum achievable softening of the material is represented by σ p − σ exp .Among these, DRX softening is the predominant type of softening observed in AZ91D magnesium alloys, and its extent can be quantified as a percentage: Materials 2024, 17, x FOR PEER REVIEW 11 of 19 type of softening observed in AZ91D magnesium alloys, and its extent can be quantified as a percentage: Jonas highlighted the determination of σs and σss values from the θ-σ curve.Using 648 K/0.001s −1 as a reference, a tangent line is applied to the θ-σ curve with σc as the reference point.The intersection of this tangent line with θ = 0 yields the value of σs, while the quadratic intersection of the θ-σ curve with θ = 0 provides the value of σss (Figure 9b).Considering various Xdrx values and utilizing Equation ( 25), σexp can be determined, with corresponding strains obtained from Figure 2 where ε = ε0.5 at Xdrx = 50%.It is evident that the true stress-strain curve of the AZ91D magnesium alloy adheres to a distinct DRX pattern, with negligible influence.Combining with Equation ( 24  According to Equation (24), one can conclude that the kd-value and βd-value is 2.38672 and 1.88619, so one can obtain the DRX model: reference point.The intersection of this tangent line with θ = 0 yields the value of σ s , while the quadratic intersection of the θ-σ curve with θ = 0 provides the value of σ ss (Figure 9b).Considering various X drx values and utilizing Equation ( 25), σ exp can be determined, with corresponding strains obtained from Figure 2 where ε = ε 0.5 at X drx = 50%.It is evident that the true stress-strain curve of the AZ91D magnesium alloy adheres to a distinct DRX pattern, with negligible influence.Combining with Equation ( 24), a linear fitting curve of ln[−ln(1 − X drx )] − ln[(ε − ε c )/ε 0.5 ] is plotted in Figure 9c, illustrating an excellent linear fitting effect.According to Equation (24), one can conclude that the k d -value and β d -value is 2.38672 and 1.88619, so one can obtain the DRX model: According to the foregoing, one can obtain Equation ( 27): From Equation ( 27), lnε 0.5 is linearly related to ln To explore the correlation between DRX volume fraction and strain, this study examined the relationship between X drx and ε, as defined in Equation (23). Figure 10 illustrates that the change in X drx can be categorized into three distinct stages as ε increases.At first, the rise in X drx is gradual, making it increasingly challenging to initiate dislocation slip, requiring higher external forces for DRX.Subsequently, X drx exhibits a rapid increase with rising ε, emphasizing the significant role of stress in driving DRX.Finally, the rate of X drx increase slows down, reaching its maximum value.Importantly, the stress needed to initiate DRX and reach this maximum value decreases with elevated temperature.This temperature-induced effect on additional dislocation generation plays a pivotal role in reducing the critical stress.Higher temperatures facilitate DRX occurrence in AZ91D magnesium alloy at lower critical values, making DRX more likely to happen.
According to Gourdet and Montheillet [38], the DRX grain size is mainly influenced by the original grain size, amount of deformation, strain rate, and alloy's deformation activation energy.As a result, the DRX average grain size can be modeled: where ε, lnd drx and 1000/T are linearly related.From the plots in Figure 11, one can find m 8 = −0.18248,Q 8 = −50,940.48J/mol, and a 8 = 713.21292.
Bringing the derived conclusion into Equation ( 29), one can derive the DRXed grain evolution model of AZ91D magnesium alloy as: requiring higher external forces for DRX.Subsequently, Xdrx exhibits a rapid increase wit rising ε, emphasizing the significant role of stress in driving DRX.Finally, the rate of Xd increase slows down, reaching its maximum value.Importantly, the stress needed to ini tiate DRX and reach this maximum value decreases with elevated temperature.This tem perature-induced effect on additional dislocation generation plays a pivotal role in reduc ing the critical stress.Higher temperatures facilitate DRX occurrence in AZ91D magne sium alloy at lower critical values, making DRX more likely to happen.According to Gourdet and Montheillet [38], the DRX grain size is mainly influence by the original grain size, amount of deformation, strain rate, and alloy's deformation ac tivation energy.As a result, the DRX average grain size can be modeled: ( ) where ddrx is the average grain size, d0 is the original grain size, Q8 is the deformation acti vation energy, a8, h8, n8, m8, and c8 is the material correlation coefficient.Since d0 is constant, a 8 d0 h8 ε n8 can be viewed as a constant a8, and taking logarithms o Equation ( 29) yields:

Grain Size under Different Deformation Conditions
Figure 12 depicts the microstructure of AZ91D magnesium alloy under a strain rate of 0.001 s −1 .Notably, an increase in temperature corresponds to a rise in the average grain size.Focusing on the highlighted area within the figure, it becomes evident that the initiation of DRX in AZ91D magnesium alloy primarily takes place at the grain boundaries.The heightened deformation temperature serves as the driving force for DRX, subsequently triggering dislocation slip and twinning.Moreover, it facilitates nucleation and migration at the grain boundaries, collectively fostering the initiation and growth of DRX and leading to an enlargement of the average grain size.Additionally, In Figure 12b,c, the white circles indicate newly produced recrystallized grains, and the arrows in the square enlargement indicate the production of twins, which play a pivotal role in the deformation of magnesium alloy, aiding in coordinating the process when dislocations impose limitations.
The heightened deformation temperature serves as the driving force for DRX, sub quently triggering dislocation slip and twinning.Moreover, it facilitates nucleation a migration at the grain boundaries, collectively fostering the initiation and growth of D and leading to an enlargement of the average grain size.Additionally, In Figure 12b,c, white circles indicate newly produced recrystallized grains, and the arrows in the squ enlargement indicate the production of twins, which play a pivotal role in the deformat of magnesium alloy, aiding in coordinating the process when dislocations impose lim tions.Figure 13 illustrates the microstructure of AZ91D magnesium alloy under varying strain rates at 673 K.A clear trend emerges from the comparison of these figures: the average grain size diminishes with an increase in the strain rate.This outcome can be attributed to the heightened deformation rate, which reduces the available time for the deformation process.Consequently, DRX lacks sufficient time to exert its softening effect, and dislocation plugging at the grain boundaries is unable to facilitate adequate grain growth.As the strain rate escalates, the concentration of dislocation plugging stress intensifies, leading to elevated nucleation rates and an increased number of new nuclei.At excessively high strain rates, the distortion energy within the recrystallized grain rises.The dislocation density within the grain reaches the threshold for initiating a new round of recrystallization before the ongoing process completes, resulting in incomplete recrystallization, as evident in Figure 13d.The dashed line in the square localized magnification shows the new grains created around the original grains, and the double-arrowed straight line indicates that the grains are compressed in the sub-direction.In this incomplete recrystallization scenario, the original grains become flattened along the direction of hot compression, and numerous new nuclei form at the grain boundaries.Unfortunately, these new nuclei lack sufficient time to mature, as highlighted in the white box, illustrating a distinctive phenomenon arising from the constraints imposed by the rapid strain rate.
shows the new grains created around the original grains, and the double-arrowed strai line indicates that the grains are compressed in the sub-direction.In this incomplete crystallization scenario, the original grains become flattened along the direction of compression, and numerous new nuclei form at the grain boundaries.Unfortunate these new nuclei lack sufficient time to mature, as highlighted in the white box, illustrat a distinctive phenomenon arising from the constraints imposed by the rapid strain rat

Finite Element Response of a DRX Model with Experimental Verification
To validate the accuracy of the previously derived hot deformation constitutive equation and DRX model for AZ91D magnesium alloy, the acquired data were utilized in Deform-3D to analyze the evolution of grain size during the hot compression process.Figure 14a,b present simulated microstructure charts depicting grain evolution under 0.1 s −1 at 648 K and 673 K, respectively.It is crucial to note that this evolution reflects the general characteristics of the grains and does not signify changes in the original grains.The different colors are intended to distinguish the grains and do not represent grain orientation.A comparative analysis of these figures reveals that with an increase in the deformation temperature, the grain size also enlarges, in accordance with the fundamental law governing grain evolution.Moreover, their grain sizes are approximately 4.78 µm and 9.65 µm, which closely align with the grain sizes listed in Table 3 (5.62 µm and 8.69 µm).Table 3 presents the average grain size of the alloy as predicted by Deform-3D simulation (marked S) and as observed in the metallographic experiment (marked E).A comparison of these values indicates that the difference between the experimental and predicted results is minimal.This demonstrates that the DRX model employed in this study is highly accurate and capable of effectively predicting the DRX behavior.To explore the evolution process of recrystallized grains, TEM canonical analysis was performed on the grain boundaries of the alloy's recrystallized grains.As illustrated in Figure 15c, throughout the deformation process, the alloy's coarse grains initiate fracture owing to deformation.Localized regions of grain boundaries accumulate a significant number of dislocations, denoted by arrows, due to large strain energy.In this energetically unstable state, certain dislocations undergo rearrangement, forming new dislocation networks and walls proximal to grain boundaries through spontaneous interactions.Concurrently, some dislocation walls gradually transform into subgrain boundaries via migration and merging, giving rise to subgrains.Owing to the high energy and instability of grain boundaries in subgrains, as deformation progresses, there is a propensity for subgrains to merge and rotate, resulting in the formation of new undistorted grains, termed DRX grains.Figure 15d presents an enlarged microstructure image of the region captured under deformation conditions in Figure 15a.Notably, small grains proximate to the larger ones, as indicated, can be observed.In Figure 15d, marked as A-grain, are small DRXed grains forming within areas of pronounced aggregation.This observation suggests a correlation between the nucleation of DRXed grains and the continuous absorption of dislocations along these subgrain boundaries, categorizing the mechanism as CDRX.Correspondingly, in the figure, distinctly labeled as B-grain and C-grain, are typical DRXed grains that are entirely isolated and enclosed by subgrain boundaries, further substantiating that an in situ transformation from subgrains to DRXed grains can be achieved by the continuous increment of large-angle grain boundaries.
grains and the continuous absorption of dislocations along these subgrain boundaries, categorizing the mechanism as CDRX.Correspondingly, in the figure, distinctly labeled as B-grain and C-grain, are typical DRXed grains that are entirely isolated and enclosed by subgrain boundaries, further substantiating that an in situ transformation from subgrains to DRXed grains can be achieved by the continuous increment of large-angle grain boundaries.

Conclusions
In this study, a comprehensive analysis of the flow stress characteristics and microstructure evolution of semi-solid forged AZ91D magnesium alloy is conducted through

Conclusions
In this study, a comprehensive analysis of the flow stress characteristics and microstructure evolution of semi-solid forged AZ91D magnesium alloy is conducted through experimental and simulation technique.This approach serves to validate the model's accuracy and visualize the changing parameters.The integrated system, encompassing physical experimentation, an ontological model, DRX behavior, and finite element simulation, underscores the significance of optimizing the forming process, controlling grain size, and enhancing overall performance.The key conclusions derived from the investigation are outlined as follows: (1) The semi-solid forged AZ91D magnesium exhibits a DRX softening mechanism during high-temperature deformation.The rheological stress of the alloy decreases with increasing deformation temperature and decreasing strain rate.And activation energy is calculated to be 156.814kJ/mol; the constitutive equation for the alloy is:   (3) Under a constant strain rate of 0.001 s −1 , elevated temperatures corresponded to larger average grain sizes.This effect was driven by increased deformation temperatures, creating favorable conditions for DRX.The higher temperature also activated intergranular dislocation slipping and twinning, promoting DRX nucleation and growth at grain boundaries, ultimately leading to larger average grain sizes.At a consistent temperature of 400 • C, an increase in strain rates resulted in a reduction of the average grain size.Higher deformation rates limited the available time for complete DRX softening, impeding grain growth at grain boundaries.However, higher strain rates induced stress concentration due to dislocation plugging, resulting in increased nucleation rates and the formation of more new nuclei.Extremely high strain rates caused original grains to flatten along the direction of hot compression, generating numerous new nuclei at grain boundaries, albeit with insufficient time for their growth.
(4) The average grain size, obtained through calculated equations, was integrated into Deform-3D to construct a material library for simulating the hot compression process.This simulation-based grain size was then compared with experimental results, thereby confirming the precision and accuracy of the calculated constitutive equations and DRX model for AZ91D magnesium alloy.

Figure 1 .
Figure 1.Schematic diagram of casting AZ91D magnesium alloy (a) and its metallography (b); schematic diagram of compressed alloys (c) and its macro-structure (d).

Figure 1 .
Figure 1.Schematic diagram of casting AZ91D magnesium alloy (a) and its metallography (b); schematic diagram of compressed alloys (c) and its macro-structure (d).

Figure 6
Figure 6 depicts a notable decrease in θ as σ values rise, marked by a distinct inflection point in the curve, indicating the initiation of the DRX process.Concurrently, the corresponding σ-value represents the critical condition for DRX.

Figure 6
Figure 6 depicts a notable decrease in θ as σ values rise, marked by a distinct inflection point in the curve, indicating the initiation of the DRX process.Concurrently, the corresponding σ-value represents the critical condition for DRX.Materials 2024, 17, x FOR PEER REVIEW 9 of 19

Figure 8 .
Figure 8. Linear fit between critical stress and strain and peak stress and strain (a) σ c − σ p , (b) ε c − ε p .

Figure 15
Figure 15 displays the microstructure for the compressed alloy with 673 K/1 s −1 .It is evident from this graph that the simulation results closely align with the actual observations, as shown in Figure 15a,b.Statistical analysis reveals that the grain size is approximately 4.89 µm, a value very close to the experimentally determined (4.51 µm) obtained from metallographic experiments.This close correspondence underscores the high accuracy of the model employed in this study.To explore the evolution process of recrystallized grains, TEM canonical analysis was performed on the grain boundaries of the alloy's

Figure 15
Figure15displays the microstructure for the compressed alloy with 673 K/1 s −1 .It is evident from this graph that the simulation results closely align with the actual observations, as shown in Figure15a,b.Statistical analysis reveals that the grain size is approximately 4.89 µm, a value very close to the experimentally determined (4.51 µm) obtained from metallographic experiments.This close correspondence underscores the high accuracy of the model employed in this study.To explore the evolution process of recrystallized grains, TEM canonical analysis was performed on the grain boundaries of the alloy's recrystallized grains.As illustrated in Figure15c, throughout the deformation process, the alloy's coarse grains initiate fracture owing to deformation.Localized regions of grain boundaries accumulate a significant number of dislocations, denoted by arrows, due to large strain energy.In this energetically unstable state, certain dislocations undergo rearrangement, forming new dislocation networks and walls proximal to grain boundaries through spontaneous interactions.Concurrently, some dislocation walls gradually transform into subgrain boundaries via migration and merging, giving rise to subgrains.Owing to the high energy and instability of grain boundaries in subgrains, as deformation progresses, there is a propensity for subgrains to merge and rotate, resulting in the formation of new undistorted grains, termed DRX grains.Figure15dpresents an enlarged microstructure image of the region captured under deformation conditions in Figure15a.Notably, small grains proximate to the larger ones, as indicated, can be observed.In Figure15d, marked as A-grain, are small DRXed grains forming within areas of pronounced aggregation.This observation suggests a correlation between the nucleation of DRXed grains and the continuous absorption of dislocations along these subgrain boundaries, categorizing the mechanism as CDRX.Correspondingly, in the figure, distinctly labeled as B-grain and C-grain, are typical DRXed grains that are entirely isolated and enclosed by subgrain boundaries, further substantiating that an in situ transformation from subgrains to DRXed grains can be achieved by the continuous increment of large-angle grain boundaries.

( 2 )
When the semi-solid forged AZ91D alloy reaches the critical strain (ε c ), dislocation accumulation reaches a critical value, leading to the occurrence of DRX.The critical strain model reflects the difficulty level of DRX.The DRX volume fraction indicates the progress of the alloy's DRX process and determines the ultimate quality of the alloy's mechanical properties.The grain size model reflects the growth trend of DRX grains.The DRX model of the semi-solid forged AZ91D magnesium alloy is calculated.Critical strain model:ε p = 0.004569 .ε 0.2662 exp(20514.87814/RT)σ c = 0.92915σ p ε c = 0.66208ε pVolume fraction model:X drx = 1 − exp −1.88619 ε−0.66208ε p ε 0.5 2.38672 ε 0.5 = 4.94040 × 10 −4 .ε 0.09446 exp(34857.35954/RT)Grain size model: d drx = 713.21292 .ε −0.18248 exp(−50940.48/RT)

Table 2 . Physical parameters of AZ91D magnesium alloy in Deform-3D [12]. Parameters Young's Modulus Poisson's Ratio Thermal Diffusion Coefficient Radiation Intensity Temperature Parameters
33208 × 10 10 (9) exp(156814/RT)(9) drx is the average grain size, d 0 is the original grain size, Q 8 is the deformation activation energy, a 8 , h 8 , n 8 , m 8 , and c 8 is the material correlation coefficient.Since d 0 is constant, a 8 d 0 h8 ε n8 can be viewed as a constant a 8 , and taking logarithms of Equation (29) yields: ln d drx = ln a 8 + m 8 ln